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O ■ ABSTRACT 

The distortion field defined by the ellipticities of galaxy shapes as projected on the sky can be uniquely 
decomposed into a gradient and a curl component. If the observed ellipticities are induced by weak 
gravitational lensing, then the distortion field is curl free. Here we show that, in contrast, the distortion 
field resulting from intrinsic spin alignments is not curl free. This provides a powerful discriminant 
between lensing and intrinsic contributions to observed ellipticity correlations. We also show how these 
contributions can be disentangled statistically from the ellipticity correlations or computed locally from 
jy-j . circular integrals of the ellipticity field. This allows for an unambiguous detection of intrinsic galaxy 

alignments in the data. When the distortions are dominated by lensing, as occurs at high redshifts, the 
decomposition provides a valuable tool for understanding properties of the noise and systematic errors. 
These techniques can be applied equally well to the polarization of the microwave background, where it 
can be used to separate curl-free scalar perturbations from those produced by gravity waves or defects. 

CO , 1. INTRODUCTION 

The shapes of distant galaxies can be distorted due to gravitational lensing by the intervening matter distribution 
■ (Gunn 1967; see e.g. Bartelmann & Schneider 1999 for a recent review and further references). While the distortions 
| are usually small, a signal can still be detected statistically since neighboring objects will be deformed in a similar way, 
. thereby producing measurable correlations in galaxy shapes. Recently, shape distortions of order 1% have been detected 
by several groups in deep galaxy surveys on scales from one up to several arc minutes (van Waerbeke et al. 2000, Wittman 
et al. 2000, Bacon, Refregier & Ellis 2000, Kaiser, Wilson & Luppino 2000). The amplitude of these distortions appears 
consistent with that predicted from weak lensing by large scale structure. 

Weak lensing is not the only possible source of galaxy shape correlations; these can also arise between physically close 
galaxies as a consequence of the galaxy formation process or subsequent interactions. One possible mechanism for this 
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' is the coupling of galaxy spins: galaxy disks tend to be oriented perpendicular to their angular momentum vectors, so 
' angular momentum couplings of neighbors will induce alignments in the projected galaxy shapes. This will also be true 
£> . but to a lesser extent for elliptical galaxies if they rotate along their shortest axis. The amplitude of the expected shape 
' correlations from angular momentum couplings have recently been studied in numerical simulations by Heavens et al. 
rS | (2000) and analytically by Crittenden et al. (CNPT 2000), while alternative mechanisms for intrinsic correlations have 
j^j ■ also been suggested (Croft and Metzler 2000; Catelan et al. 2000). Intrinsic correlations will be especially important in 
relatively shallow surveys such as the 2dF or the Sloan Digital Sky survey, where the median redshift is <C 1. Evidence 
for the existence of such intrinsic correlations in the nearby universe has been recently presented by Pen, Lee & Seljak 
(2000) and Brown et al. (2000), looking in the Tully and the Super COSMOS galaxy catalogs, respectively. 

One way of disentangling the intrinsic shape correlations from those induced by weak lensing is to examine the patterns 
of the average galaxy shapes. The shapes of galaxies are typically described by two degrees of freedom; their average 
ellipticity and orientation. The distribution of galaxy shapes can be described by a symmetric and traceless 2D tensor 
field. In general, any such tensor field can be written as a sum of two terms, one of which is curl-free and the other is 
divergence-free. In analogy with the radiation field in electromagnetism, these are usually referred to as the 'electric' (E) 
and the 'magnetic' (B) component, respectively. 

Lensing by a point mass will create a tangential, curl-free distortion pattern. The most general distortion field produced 
by lensing will be a linear superposition of such patterns and will also be a curl-free, (i.e. -E-type) field (Kaiser 1992, 
Stebbins 1996). However, as we show below, the distortion field resulting from intrinsic spin alignments has E and B- 
type modes of the same order of magnitude. This property will enable us to uniquely disentangle angular momentum 
correlations and to subtract their contribution from a measured distortion field in order to more accurately compute and 
isolate the distortions induced by lensing alone. When the deformations are dominated by lensing, the E~B decomposition 
can improve the signal to noise level, since noise and other systematic effects are expected to contribute to both the E 
and B channels. Current surveys have measured the sum of E and B powers, thus doubling the noise power relative to 
the decomposition strategy that we propose here. 

The expected amplitude of galaxy shape correlations is fairly small. In order to measure it given the large scatter in the 
intrinsic shapes of galaxies, many galaxies must be observed. At present, observations of mean ellipticities are dominated 
by the intrinsic scatter. Direct decomposition into E and B modes is a non-local operation, requiring derivatives of these 
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noisy observations and so is very problematic. Here we show how the correlation functions of the observable ellipticities 
can be directly converted into correlations of the E and -B-modes. 

It is also useful to have locally defined quantities which reflect the E and B decomposition. Kaiser et al. (1994) and 
Schneider et al. (1998) looked at this issue in the context of lensing and developed a statistic known as the 'aperture 
mass' which enabled them to put a lower bound on the projected mass in a localized region of the sky. More generally 
this statistic gives a direct, local measure of the electric contribution to the distortion field, and a similar observable can 
be evaluated to measure the magnetic component. We develop this formalism here and relate the correlations of these 
local E and B estimators to those of the observable ellipticities. 

This paper is organized as follows. We begin by defining the E — B decomposition in terms of observed ellipticities, and 
demonstrate that the lensing distortions are curl free, whereas those due to intrinsic alignments are not. In Section 3, we 
discuss estimators of the E — B correlations and their relation to correlators of the ellipticities. In Section 4, we define 
the local E and B measures and calculate their correlations. We conclude in sections 5-6 with a cookbook style summary 
on how the E — B decomposition can be derived from the statistical weak lensing surveys. 

2. DECOMPOSITION OF THE DISTORTION FIELD 

The projected shape of a galaxy on the sky can be approximated by an ellipse with semi-axes a and b (a > b), of which 
the major axis makes an angle tp with respect to the ir-axis of the chosen coordinate system. It can then be concisely 
written as a complex number, 

where e + — |e|cos(2?/>) and e x = |e|sin(2^). The phase dependence oc e 2 ' 1 ^ expresses the fact that the ellipticity is 
invariant under a rotation over n radians. Note that the two components e x and e + are analogous to the Q and U Stokes 
parameters for linearly polarized light. Given a distribution of galaxies on the sky with measured ellipticities, the complex 
scalar ellipticity field defines a traceless, symmetric 2x2 tensor field, 



[l]ab = 



e x -e_i 



(2) 



We have approximated the sky as flat and followed the derivation of Kamionkowski et al. (1998). See Stebbins (1996) for 
the generalization to a curved sky. 

The shear field Eq. (|J) can be written in terms of a gradient or 'E' -piece and a curl or pseudo-scalar '£?' -piece (Stebbins 
1996) by introducing two scalar functions $e and <&b, 

7ab (x) = (d a d b - i(UV 2 ) * B (x) + ^(e cb d a d c + e ca d c d b ) $ B (x), (3) 

where e ab is the anti-symmetric tensor. Each component of the ellipticity field can be written as a function of $e and 
$b as 

e+ = Ixx = -lyy = -jipA - dydy) $ B (x) - d x d y <£> B (x) 

£x = lyx = Ixy = d x d y $ E (x) + ~(d x d x - dydy) <J> s (x). (4) 

The E and B parts can be extracted explicitly from the shear tensor by applying the V operator, 

V A <S> E = 2d a d blab ; V 4 <5> B = 2e ab d a d clbc . (5) 

The relation between the functions and <&b and the projected gravitational potential will become obvious in what 
follows. It is useful to perform the E — B decomposition in terms of variables that have the same dimension as the 
measured ellipticities (Kamionkowski et al. 1998): 7# = iV 2 <I>£ and jb = ^V 2 <&b- These are related to the ellipticities 
by 

V 2 7 B = d a d b j ab = (d x d x - d y d y )e + + 2d x d y e x 
V 2 7s = e a b d a d c j bc = (d x d x - d y d y )e x - 2d x d y e + . (6) 

Since only the derivatives of and 73 are defined in terms of the ellipticity field, 7b (x) and 75 (x) are ambiguous up to a 
constant and linear gradient term. In the same way, constant and linear gradient terms in the measured ellipticities should 
not impact the E — B decomposition. To see this, consider an ellipticity field where e+(x) = x and £x( x ) = 0. This can 
either be a consequence of a pure E'-mode ($£ = x 3 /3; $b = 0), the result of a pure B-mode ($£ = 0; <f>s = x 2 y/2+y 3 /6), 
or a linear combination of the two. 

A rotation of the basis axes translates e+ into e x and vice versa, but does not affect the E — B decomposition. In 
particular, the ellipticity measured in a basis which is at an angle ip relative to the original basis is given by: 

e' + = e + cos 2ip — e x sin 2ip ; e x = e + sin 2ip + e x cos 2(p. (7) 

Thus, a global rotation of 7r/4 transforms e' + — — e x ; e' x = e + , but since the position vectors are also rotated, the E — B 
decomposition remains invariant. However, one can also take an ellipticity field and rotate each ellipticity individually by 
7r/4, keeping the position vectors fixed. This new ellipticity map has the E and B modes of the original map interchanged: 
Ye = -lB;i B = IE- 
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2.1. Distortions due to lensing 

In the case of gravitational lensing, the resultant distortion field 7 can be written in terms of a gravitational deflection 
potential ip as (e.g. Bartelmann and Schneider 1999) 

7ab(x) = (d a d b - ^ afc V 2 )^(x) . (8) 

The deflection potential tjj is a convolution over the projected surface mass density k(x), il>(x) = — J dx'rt(x') In |x — x'|. 
Comparing this expression to the E — B decomposition of Eq. (p|) we can identify ^(x) = ip(x) and = 0. Thus for 
the shear field induced by lensing, the E-mode is related to K, the projected surface mass density in units of the critical 
surface mass density for a given configuration of source and lens, and the B-mode is identically zero (as was discussed by 
Kaiser 1995; Kamionkowski et al. 1998). Note that weak lensing only approximately gives pure E-modes, as B-modes 
may arise when the light is bent in more than one scattering event. However, these B-modes arise at higher order and 
are suppressed relative to the E-modes. 

If one measures the E-mode contribution of a given map and then rotates every measured ellipticity by 7r/4 and repeats 
the same measurement, one obtains an estimate of the B-contribution which should be consistent with zero for a pure 
lensing signal. Therefore, the absence of B-modes naturally provides a robust test for isolating the lensing component of 
the distortion field and provides an estimate of the noise level of the data (Kaiser 1992). 

Note that lensing is not the only possible source of curl-free correlations. If the shapes of galaxies are primarily 
determined by tidal stretching, this would lead to intrinsic correlations (Catelan, Kamionkowski & Blandford, 2000; Croft 
& Metzler, 2000). For these shape distortions, the observed ellipticities are also linear in the tidal field, leading to pure 
electric modes just as in lensing. For spiral galaxies which have had many dynamical times to evolve, tidal stretching is 
likely to be small compared to the contribution from spin alignments. Even for elliptical galaxies, a bulk rotation of as 
small as 1 km/sec would erase the galaxy's original alignment. Almost all observed ellipticals rotate faster than that, so 
the shape-shear alignment in Catelan et al. (2000) is unlikely to be observable. However, it could be significant for larger 
objects like clusters which are dynamically much younger. 

2.2. Distortions due to angular momenta alignments 

Shape correlations between galaxies can also arise from alignments in the direction of their angular momenta. While this 
is particularly true for spiral galaxies, given the assumption that their disks are perpendicular to the angular momentum 
vectors, it is true to a lesser extent for elliptical galaxies as well. Ellipticals probably rotate about their intermediate axis 
(Dubinski 1992), so averaging over all statistical randomized alignments, the average major axis is perpendicular to the 
angular momentum vector, just like a spiral galaxy. The strength of the correlation signal has recently been studied in 
numerical simulations by Heavens et al. (2000). We have recently attempted to model these theoretically (CNPT 2000) 
by assuming angular momentum is induced by tidal torques. Following the formalism developed by Catelan and Theuns 
(1996), the correlations can be calculated for Gaussian initial fluctuations using linear theory coupled with the Zeldovich 
approximation. 

The induced intrinsic correlations of ellipticities will primarily result from correlations in the direction of the angular 
momenta. As discussed in CNPT this implies that the ellipticities are effectively quadratic in the angular momenta, 

^+ ~(T x iTi x T^yiDiy) , £x C>C ~^Dx{D%y . (9) 



Here, Ty cx didj(f> is the shear of the gravitational potential, and T denotes the shear tensor normalized by ^/TijTij. Note 
that i and j run over three coordinates, x,y,z in contrast to above where two dimensional (projected) quantities were 
considered. The quadratic dependence on the shear in Eq. (|^) is fundamentally different from the linear one appropriate 
for the lensing case, Eq. (^), and as a result the -B-modes are non-zero. 
In particular, it is straight forward to show that 

V 7-E Txi^xxTix -\- Tyi yyTiy -\- T x i : xy^iy ~\~ Tyi^xyT-'ix ~\~ {^xi. x Tyi^y) . (10) 

and similarly 

V ^ Txi^xxTiy Tyi yyTix T^xi^xy^ix ~i~ T^yi^xy^iy- (H) 

The amplitude of the -B-modes is comparable to that of the E-modes, as is shown in Figure 1. The presence of B- 
modes provides a mechanism for disentangling the correlations which arise from lensing from those resulting from intrinsic 
alignments due to angular momentum couplings. 

This is not the sole means of disentangling intrinsic correlations from weak lensing. The observed ellipticity correlations 
from intrinsic alignments are strongest at low redshifts, while the lensing signal is larger when the sources are at higher 
redshifts. In addition, morphology distinctions will be useful, since intrinsic correlations of elliptical galaxies are smaller 
than those of spirals because they are intrinsically more round. Further possible methods for distinguishing lensing from 
intrinsic correlations are discussed in Catelan, Kamionkowski and Blandford (2000) and CNPT (2000). 

3. CORRELATION ESTIMATORS 

The decomposition into curl and gradient contributions is most straight forwardly performed in Fourier space (Kamion- 
kowski et al. 1998). However, often it is useful to consider this in real space, as many issues which might complicate 
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matters in Fourier space, such as finite field size or patchy sampling, are more easily handled in real space. In this 
section we deal with performing the decomposition statistically in real space using the ellipticity two point functions. 
This approach is particularly relevant when the observations are noise dominated, such as is the case when there are 
relatively few galaxies with which to measure the mean ellipticity, which is the case for the current surveys on small scales 
where most of the lensing signal lies. In the next section, we will address the issue of a local decomposition. 

3.1. Correlations in je and jb 

Here we will relate the correlations of the electric and magnetic shear directly to correlations of the ellipticity. The real 
space correlation function of je is related to correlation of its corresponding potential, <5>e, by 

£e = (7b(x)7b(x + r)> = i(V 2 <Mx)V 2 $ B (x + r)) = \v 4 Z E (r), (12) 

where Es(V) = (&e (x) &e (x + r)). An analogous relation holds for 75 correlations, while the cross correlation, (7_e7b)j 
is zero if the field is invariant under parity transformations. 

These correlations can be directly computed from the observed correlations in e + and e x , defined through 

C x (r,<p) = <e + (x)e + (x + r)); C 2 (r, <p) = (e x (x)e x (x + r)) (13) 

where the ensemble average is over pairs with separation r, for which the separation vector r makes an angle ip with 
respect to the chosen basis. The sum of these correlations is rotationally invariant, but their difference depends explicitly 
on the choice of orientation of the coordinate axes (Kamionkowski et al. 1998). 

The required relation between the observed correlations C\ and C2, and the E — B-correlations £e and £b follows from 
their definitions using Eq. (||) written in terms of derivatives with respect to the separation r: 

Ci{r, ip) = iv 4 [S B (r) + 3 fl (r)] + ^x[S B (r) - H B (r)] cos4<^> 

C a (r,<p) = ^V 4 [S B (r) +S fl (r)] - ^x[S B (r) - E B (r)] cos4^, (14) 

where V 4 = 8D 2 + 8r 2 D 3 + r 4 D 4 , the operator x — r 4 D 4 and D = We wish to invert these equations to find 

expressions for £e and £b in terms of the observable correlation functions. 

This inversion can be done most easily in terms of basis independent correlation functions which we shall denote £+ 
and £ x . Physically, corresponds to the correlation function (e + (x)e + (x + r)) computed in such a way that, for each 
pair of galaxies, one coordinate axis is always taken parallel to the separation vector r. A similar definition holds for £ x 
in terms of (e x e x ), while by isotropy, the expectation of the cross correlation is zero. These correlations are related by a 
rotation to C\(r, tp) and Ci(r, ip) and satisfy the relations, 

+ Cx (r) = d(r, V ) + C 2 (r, <p) ; g + (r) - [r)\ cos(4^) = d(r, <p) - C 2 (r, <p) . (15) 

In terms of these new observables, Eq. ( |l4| ) simplifies to 

4K+(r) + £x (r)] - V 4 [S £ (r) + S B (r)]; 4[€+(r) - Cx (r)] = x|3 B (r) - S B (r)]. (16) 

Initial measurements of the lensing signal have focused primarily on the variance of the magnitude of the ellipticity 
averaged over regions of a given size, and its fall off as the size of the regions is increased. The variance is simply the value 
of the correlation C\ + C2 , convolved with the appropriate window function, at zero lag. The variance has the advantage 
that it is a local quantity and is straight forward to measure. However, the measurements of the variance at different 
scales have strongly correlated errors. In addition, since C\ + C2 — |V 4 (Ss + Hb), such measurements are unable to 
distinguish _E-modes from B-modes. Therefore, it is advantageous to investigate the full correlation functions. 

The inversion of Eq. ( |I"4| ) can now be written as, 

Ur) - ^[tf(r) +£x(r)] + ^VV^+M - £ x (r)] 

Ur) = llt + (r)+ZM - ^VV'K+W -Ur)]- (17) 
An equivalent set of equations follows by applying ^V~ 4 operator to both sides of these equations. 

x V- 4 Ur) = JxV- 4 K + M +*x(r)] + i[f+(r) -C x (r)] 

xV- 4 ^(r) = I xV - 4 K+(r) +CxM] - i[C + (r) . (18) 

This form actually proves more useful in practice since, as we show below, x^ -4 is more local than V 4 x _1 - Note that 
these expressions assume statistical isotropy, that is, (e + e x ) = in the basis where the axes are aligned with the galaxy 
separation vector. 
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3.2. Evaluation of V 4 \ 1 and %V 4 

It is quite useful to take these relationships into Fourier space, and the operators V 4 and x have particularly simple 
expressions when applied to the Bessel functions which arise in a Fourier transform. In particular, 

V 4 J (fcr) = fc 4 J (fcr); xM^r) = k 4 J 4 (kr) . (19) 

When combined with the above equations, we can relate the correlation functions directly to the E and B power spectra, 
reproducing the relations Eqs. (14) and (15) of Kamionkowski et al. (1998). Also, it is straightforward to show that 
since the power spectrum is the Fourier transform of the correlation function, £_E(r) = J kdkJo(kr)PE(k) implies that 
X^- i iE{r)= jkdkHkr)P E {k). 

The operators V 4 x _1 and x V~ 4 are also most easily evaluated in Fourier space and can be shown to take the form, 

V 4 x^s(r)= /^ J °( fcr ) / r'dr'J 4 (kr')g(r') = J r'dr'g(r')g(r,r') 
xV- 4 . 9 (r) - J ^Mkr) J r> dr> J (kr')g(r') = J r> dr> g(r')G(r' ,r) , (20) 

where 

^-J Q (kr)J 4 (kr'). (21) 

Note that since the k integral is from zero to infinity, rr'Q(r,r') depends only on the ratio r/r'. This operator is not 
simply a convolution as would be the case if Q(r, r') were a function only of the difference |r — r'|. In fact, the Fourier 
space operator takes nearly the same form as the real space operator. That is, if the Fourier transform of g(r) is g(k), 
then the transform of J r'dr'g(r')Q(r,r') is J k'dk'g(k')G(k' ,k). 
Alternatively, these operators can be written in integral form, 



4 5 (r) = g(r) +4 / dr> ^ Ur 2 H ' dr> ^ 



vV 

xV- 4 5 (r) = g(r) + ^ f dr'g{r')r' - ™ f dr' g(r')r' 3 . (22) 
r Jo r Jo 

The second form is the more useful in practice, being more local. To evaluate it at some radius R, one only needs to know 
the form of g(r) for r < R, whereas the V 4 x' 1 requires knowing g(r) for r > R. 

The particular form of the first of the integral operators follows from choosing integration constants such that V 4 (S# (r) + 
Ss(r)) does not diverge at large separations. To obtain the second form, we demand that x(S^(r) — H^(r)) and its deriva- 
tives are well-behaved as r — > 0. The origin of these integration constants is the constant and linear gradient ambiguities 
in the definitions of -f E and -f B in Section 2. 

3.3. Power law solutions 

As an illustration, consider the case where the correlation functions both behave as power laws with the same index, 
£+(r) = Ar n and £ x (r) = Br n . Note that V 4 x _1 r n = f{n)r n , where f(n) = (n 2 + 6n + 8)/(n 2 - 2n). It then follows 
that 

He,b - f(n) X y-^E,B = \[{A + B)±{A- B)f(n)]r n . (23) 

Note that for certain power laws, n — — 2 or n = —4, the operator V 4 x~ 1 r n — 0, while it diverges when n = or n = 2. 
The opposite is true for the inverse operator x^~ 4 , which diverges when n = — 2 or n = — 4 and is zero for n = or 
n = 2. 

In the case when the B-modes are exactly zero, as occurs for gravitational lensing, one obtains (A + B) = (A — B)f(n). 
It follows that the ratio of the correlation functions is given by 

t it n 2 + 2n + 4 
4(n + 1) 

This expression diverges when n = —1 which therefore implies that £ x = in this case. Kaiser (1992) also considered 
power law spectra in the lensing case and presented results for a number of spectral indices. Our results agree qualitatively, 
but the agreement is not exact. Kaiser calculated these numerically and this may be the source of the discrepancies. 

3.4. Application to spin correlations 

We can apply this technique to the correlations arising from intrinsic spin couplings. The ellipticity correlations were 
calculated for the model described in CNPT (2000), and we will not go into further detail here. Using the CNPT results 
and the above expressions, we can calculate the E and B correlation functions and these are presented in Figure 1. In 
contrast to the gravitational lensing cases, the E and B modes are seen to be of comparable magnitude in this model. 

Various parameter choices have been made in the examples we have shown, but we believe the E—B decomposition to be 
largely independent of most of these. For simplicity, the figures assume the galaxies are effectively perfect disks (a = 1 in 
the notation of CNPT). Using more realistic galaxy shapes will only suppress the overall amplitude of the correlations. We 
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Fig. 1. — The E and B-mode correlation functions for intrinsic spin correlations in the model of CNPT (2000). The amplitude of the 
correlations are determined by the parameters a and a, which have been taken to be unity for simplicity. The mean redshift of the sources 
was taken to be z m = 0.1 and the density correlations were taken to fall off as r~ 1 in the left figure and as r -3 / 2 in the right figure. Also 
plotted are the differences between £e and f+i which is the same as the differences between £b and £ x - In the left panel, the projected 
cllipticity correlations fall off as and £g and are very nearly the same. This is not the case in general, as can be seen in the right 
panel where the projected correlations fall as 8~ 2 . However, the 0~ 2 is also special in that t;g and £g are nearly identical. 



have also assumed that the angular momentum directions correspond with those that would be predicted by linear theory 
(a = 1.) Non- linear evolution may affect the direction of a galaxy's angular momentum, but as long as these changes are 
not coherent then they will only suppress the overall correlation amplitudes. Finally, we have assumed a mean redshift 
for the sources of z m = 0.1. Changing this will change the angular scale at which a given level of correlations are seen, 
but it should not affect the nature of the E — B decomposition. 

As discussed above, one factor which could affect the E — B decomposition is the rate at which ellipticity correlations 
drop off. The case when the density correlation falls off as r _1 implies that for large separations, the projected ellipticities 
fall off as 8~ 1 , where 9 is the angular separation. Since /(— 1) = 1, £e will be the same as £+, and £g will be the same 
as £x . This is shown in the left panel of Figure 1. Also shown in the right panel is the case when the density correlation 
falls off as r -3 / 2 and the projected ellipticity drops as 8~ 2 . In this case, /(— 2) = and £,e ~ (,b — +Cx)- Note that 
here the x modes are actually anti-correlated at large separations. 

4. LOCAL CORRELATORS 

Here we will consider local estimators of the E and B modes. These are generalizations of the 'mass aperture' formalism 
of the pure lensing case, proposed by Kaiser et al. (1994) and developed further by Schneider et al. (1998) and applied to 
CMB polarization by Seljak and Zaldarriaga (1998). They showed that a convolution of the tangential shear with a given 
wavelet provided a measure of the projected mass convolved with a related wavelet. More generally, the integrals of the 
tangential shape distortions can be directly related to the local 'electric' distortion. Thus -E-modes are associated with 
either tangential or radial patterns, as shown in Figure 2 (left panel.) B-modes are also related to the circular distortion 
pattern, but these have an associated 'handedness' or orientation, as shown in the right panel of Figure 2. 

The local correlators are most easily defined by considering polar coordinates about a given point. The E-mo&e is 
related to the tangential shear, which is the local e + (or Q Stokes parameter) in a coordinate system defined by the radial 
vector to that point. The B-mode corresponds analogously with the local e x mode (or U Stokes parameter) in the same 
basis and can be thought of as a ir/4 shear. These quantities are related to the original cllipticity field by a ip dependent 
rotation: 

-ft — e+ cos(2<^) + e x sin(2<^); jjl = e x cos(2tp) — e + sin(2<^). (25) 

Figure 2 shows modes where 7t(<£>) and 721 (ip) are independent of ip. 

It is possible to show that circular integrals of 74 and 7^ are directly related to the and 73 contained interior to the 
circle. Using the relations of Eq. (Q) in polar coordinates, one can show 

- I (&L _ I A _ J_J^ $ -(JLlJL\i> - 1 ( d% 19 1 92 (JLIJL\$ (26) 

' 2 \8t 2 r dr r 2 dip 2 ) E \dr r dip J B: ^ 4 2 \dr 2 r dr r 2 dip 2 ) B \dr r dip ) E 

We can use the polar form of V 2 to relate the derivatives of the potential to 7^. If we integrate these equations over ip, 
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B - modes 



Fig. 2. — Local representations of E (gradient) and B (curl) modes. E-modes are either tangential or radial, depending on their sign. 
B-modes can be oriented in either a clockwise or counter-clockwise (shown) direction. Lensing generally brings about only E-modes, while 
noise and angular momentum correlations can generate both. Local estimators of the E and B-modes can be found by doing a radial weighting 
of these circular integrals (Kaiser et al. 1994). 
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o 

2n 



<^7t(r) + 2 / ^r— I ^Tt(r') (27) 
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A similar relation holds when replacing j E — ► 7b and 74 — > 7^ . 

We can convolve these relations with a compensated filter (r) . This filter may be arbitrary, but we will require that 
/ d 2 rU(r) = 0. For example, one might take it to have the shape of a Mexican hat. Multiplying both sides by rU(r) and 
integrating over r one obtains the local estimators 

1 / d 2 r lE {r)U{r) = J d 2 r lt {r)Q{r) 

(28) 



T B = J d 2 r lB {T)U{r) = J d 2 r 7f (r)Q(r) 



where it follows by integrating by parts that Q(r) = U(r) — pr f Q r'U{r')dr' '. U{r) can be taken to be zero outside a 
given radius, so these relations become purely local. Thus we have a local measure of the E — B decomposition related 
solely to the ellipticity in that region. In the absence of instrumental and sampling noise, lensing predicts that Tb will 
be identically zero for any point on the sky. 

The correlations of these local measures can be computed as follows, 

d 2 k. 



(T E (0)T E (R)) = I ^U 2 {k)e 



ik-R 



d 2 r- [K+(r) + (r)) + VV 1 {£+(r) - £x (r)}] W(\r + R|) 



= r / d 2 r[^ + (7 



• £ x (r)]W(|r + R|) + i / ,Pr[U ( r) - ; x (Vi]l Vi r + R|). 



(29) 



where we used Eq. (|l7|), and defined VV(fe) = U 2 (k), so that W(r) is the convolution oiU{r) with itself. The function 
VV(r) can be shown to be a convolution of W(r) and ^(r, r'), i.e. 



7'dryV(r)V 4 x : g(r) = / rdrW(r) / r'dr'Q(r,r')g(r') = / r'dr'fl'(r / )W(r / ) 



(30) 



From this it follows that VV(r') = J rdrW(r)5(r, r') = xV 4 W(r'). The corresponding expression for the _B-mode is 
(Tb(0)T b (R)) =\l d 2 r[C + (r)+^(r)}W(\r 



Rl 



^r[C+(r)-ex(r)]W(|r + R|) 



(31) 



The variances of the E or B field smoothed with a given window W (r) is obtained by setting the separation R = 0. Note 
that as in the previous section, we are implicitly incorporating statistical isotropy in these correlation expressions. 

For concreteness, it is useful to consider a simple example of a wavelet shape. Following work by van Waerbeke (1998). 
assume the radial function to have the form of a Mexican hat wavelet, which is the derivative of a Gaussian function, 
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U(r) — o~ 2 (1 - r 2 /2a 2 ) exp(-r 2 /2cr 2 ) and its Fourier transform is simply U(k) = \k 2 a 2 e~ k2 ^ 'I 2 . For this particular 
choice, the convolution of U{r) with itself is W(r) = (l/2cr 2 ) [2 — r 2 ja 2 + r 4 /16a 4 ] exp(— r 2 /4a 2 ). Finally, using the fact 
that the Fourier representation of V -4 is fc~ 4 , we have, 

J_ ( L 

2a 2 \Aa 2 



Mr) = X 



e" r /4<J . (32) 



This particular wavelet has the advantages that it is simple, analytic and very compact, falling off exponentially at large 
distances. 

5. DATA ANALYSIS 

In this section we address the direct analysis of real survey data. We first consider the extraction of the E and B 
correlators. For this purpose, one only needs to measure the two pairwise ellipticity correlation functions (defined 
before Eqn (15)) for all pairs of galaxies as a function of separation, which has been done by van Waerbeke et al, (2000) 
and does not depend on the geometry of the survey or its boundary shapes. This correlation function (plus its error bars 
and the covariance matrix of errors) contains all the second order statistics of the map, and is a complete and optimal 
two point description. 

We now have two functions, both of which contain noise, and our goal is to apply a rotation such that one function 
contains lensing signal and the other no lensing signal, but which will give an estimate of contamination from noise 
and intrinsic correlations. The current analyses effective have added the two correlations, which adds a function which 
contains lensing to one which contains no lensing but an equal amount of noise, which doubles the amount of noise we 
have. Instead, we can define £'(r) = 2/r 2 f Q r [£ + (r') + £x (r'))r'dr' — 6/r 4 J^[^ + (r') + £ x (r')]r' 3 dr' . We can now derive pure 
£J-type and B-type correlators which depend only on correlations at separations less than r, 



xV- 4 £ £ (r)= / kdkMkr)P E (k)=t + (r)+t'(r) 

xV" 4 £ B (r) - / kdkMkr)P B (k) = £x(r) +£'(r) (33) 



In the case of pure weak lensing, the B-type correlator should be consistent with pure noise, while £g contains all the 
lensing signal, and only half the noise. We have achieved the correlation function analogy of Kaiser's 45 degree rotation: 
rotating all the images by 45 degrees swaps £g and £b. 

One can also obtain expressions for the variances of the fields smoothed by a tophat filter with radius R. This is done 
by convolving the correlation functions with a window which is proportional to the area of overlap between two circles of 
radius R and separation r, 

(7%{R))th = j\dr^(r) (2R 2 cos^ 1 (r/2i?) - r^R 2 - r 2 /4) . (34) 

This again contains half as much noise power as the standard procedure of actually convolving the map and computing 
its variance, which is what all work to date has performed. 

This decomposition was performed directly on the correlator, and there exists no transformation on the shear map such 
that the correlation function of the transformed shear map is given by Eqn. d33|). Similarly, Eqn. (|34|) is not identically 
equal to the tophat smoothed map variance. If one actually wants to make a map, the aperture shear decomposition Eqn. 
(28) allows one to make a local decomposed map, whose variance one can then measure. Note that in making a map, one 
looses information near the boundaries, and variation in source counts leads to inhomogeneous signal-to-noise, decreasing 
the overall sensitivity to measuring the true correlation function. Luckily the aperture mass approach also allows a direct 
computation of the same quantities from the correlation functions, as shown in Eqn. (31). 

6. CONCLUSIONS 

Here we have investigated the general decomposition of flat two dimensional spin-2 fields into so-called electric (gradient) 
and magnetic (curl) components. While this decomposition involves derivatives and is intrinsically non-local, we have 
shown how local correlations of the electric and magnetic components can be found given correlations in the components 
of the ellipticity (Equations 17 and 18). For the case of power law correlations, this implies a relationship between the 
spectral index of the correlations and the relative amplitude of the different types of ellipticity correlations. 

In addition, following Kaiser et al. (1994) we have shown how local estimators for the electric and magnetic modes can 
be constructed from circular integrals of the tangential and 7r/4 (or Q and U Stokes parameters) distortions respectively 
(Figure 2). We calculated correlations of these local estimators and related them to the ellipticity correlations. 

This decomposition has important consequences when applied to the projected shapes of galaxies. Gravitational lensing 
primarily produces electric modes, as does the tidal stretching of galaxies. However, as we have shown here, angular mo- 
mentum couplings produce E and B-modes in comparable amounts and one might expect that noise, telescope distortions 
and other sources of systematic errors will produce curl modes as well. 

Thus the presence of B-modes will be useful for disentangling intrinsic correlations caused by angular momentum 
couplings from those induced by cosmic shear and from gravitational lensing. Even if lensing distortions dominate, this 
decomposition will be useful as a means of estimating the levels of noise and systematic errors in the observations. In 
addition, it provides a means of reducing noise levels of lensing observations by a factor of \[2 (Kaiser 1995.) 
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The prospects for isolating the contribution of -B-modes using these local correlators are promising with several of the 
ongoing shallow redshift and imaging surveys presently taking data. These include the Sloan DSS (www.sdss.org), 2dF 
(www.ast.cam.ac.uk/AAD/2df), 2MASS (pegasus.phast.umass.edu), DEEP (dls.bell-labs.com) and the INT wide 
field survey (www. ast . cam. ac . uk/^wf csur/). The low median redshift implies minimal contamination from the lensing 
signal and therefore an improvement in the signal to noise of the extraction. 

Finally, most of these considerations apply equally well to the imminent observations of CMB polarization (e.g. 
Kamionkowski, Kosowsky & Stebbins 1997; Zaldarriaga & Seljak 1997). In this case, scalar fluctuations induce only 
-E-modes, while noise and gravitational radiation induce both E and _B-modes. As in the case of galaxy shapes, these 
observations will be initially noise dominated, so these kinds of correlation analyses will be essential. The lessons we learn 
from the lensing data now available will be directly applicable to the polarization data when they become available in a 
few years time. 



We thank Ben Metcalf, Neil Turok and Ludo van Waerbeke for useful conversations. RC and TT acknowledge PPARC 
for the award of an Advanced and a post-doctoral fellowship, respectively. 
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